第 40 卷 第 6 期 应 用 力 学 学 报 Vol. 40 No.6 
2023 年 12 月 Chinese Journal of Applied Mechanics Dec. 2023 


文章 编号 :1000- 4939 (2023)06- 1412- 09 


基于 CFD 的 两 林 式 请 次 机 夫人 
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摘要 :针对 某 两 栖 式 清 浴 机 器 人 ,基于 计算 流体 力学 (computational fluid dynamics ,CFD ) 分 析 其 静 
浮 态 和 航行 阻力 。 采 用 重 登 网 格 方法 计算 静 浮 态 下 的 吃水 深度 和 纵 倾 角 ; 分 别 采 用 静态 网 格 方法 
和 重合 网 格 方法 预测 两 栖 式 清洲 机 器 人 静水 直 航 工 况 的 航行 阻力 ,并 对 这 两 种 方法 进行 对 比 。 研 
究 结 果 表 明 : 采 用 重合 网 格 的 方法 可 以 有 效 计算 两 栖 式 清 淤 机 器 人 静 浮 态 下 的 吃水 深度 和 纵 倾 角 ; 
静态 网 格 方法 和 重 司 网 格 方法 均 能 有 效 预 测 两 栖 式 清 湾 机 器 人 水 上 航行 阻力 ,静态 网 格 方法 计算 
效率 高 ,但 仅 适 用 于 航速 较 低 的 工 况 , 而 重 胎 网 格 方法 适用 于 不 同 航速 下 的 航行 阻力 预测 ,但 计算 
效率 较 低 。 研 究 结 果 对 两 栖 式 清洲 机 器 人 的 设计 及 改进 有 较 大 参考 价值 。 
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Research on hydrodynamic characteristics of amphibious 
desilting robot based on CFD 


RUI Hongbin, ZHENG Wenzhe ,LI Lulu,CAO Wei 


(School of Mechanical and Precision Instrument Engineering ,Xian University of Technology ,710048 Xian ,China ) 


Abstract: Based on computational fluid dynamics (CFD), the floating state and sailing resistance of 
amphibious desilting robot were analyzed. The draft depth and pitch angle were obtained by calculating the 
floating state with the overlapping mesh method. The static mesh method and overlapping mesh method 
were used to predict the sailing resistance of amphibious desilting robot and the two methods were 
compared. The results show that the draft depth and pitch angle of the amphibious desilting robot can be 
calculated effectively by using the overlapping mesh method. Both the static grid method and the 
overlapping mesh method can predict the sailing resistance of amphibious desilting robot. The static mesh 
method has high efficiency, but can only apply to the Situation of the low sailing speed. While the 
overlapping mesh method is suitable for sailing resistance prediction under different sailing speed ,it has 
low efficiency. The results are of great reference value to the design and improvement of amphibious 


desilting robot. 
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近 几 年 来 ,大 规模 土地 开发 加 上 连续 暴雨 天 气 ， 
导致 我 国 河道 湖泊 的 淤泥 沉积 问题 严重 ,该 现象 大 
大 前 减 了 我 国 划 水、 排水 的 能 力 ,洪涝 灾害 时 常 发 
生 。 因 此 , 河 湖 的 政 浚 、 清 淤 工 作 得 到 高 度 的 关注 和 
重视 。 我 国 拥有 大 量 复杂 河 网 和 中 小 型 河道 , 玻 
浚 清洲 工作 复杂 且 难 度 大 ,人 工 清 浴 、 船 载 控 掘 机 
和 大 型 清 汶 船 等 不 适用 于 这 些 复杂 工 况 ,其 工作 效 
率 、 自 动 化 和 智能 化 程度 较 低 ,而 中 小 型 两 栖 式 清 汶 
机 器 人 凭借 其 较 高 的 自动 化 和 智能 化 ,在 清洲 行业 
有 广阔 的 发 展 前 景 。 

目前 两 栖 式 清洲 机 器 人 水 动力 特性 研究 较 少 ， 
考虑 到 两 栖 式 清 汶 机 如 人 作为 一 种 特殊 的 水 上 航行 
器 并 结合 其 工作 特点 ,发 现 浮 性 和 快速 性 是 其 水 动 
力 特 性 中 的 研究 重点 。 浮 性 是 两 栖 式 清洲 机 器 人 最 
基本 性 能 ,是 其 他 性 能 研究 的 基础 ,在 一 定 程 度 上 反 
映 了 机 顺 的 承载 能 力 ,其 核心 问题 是 静 浮 态 的 计算 ; 
快速 性 对 于 机 需 在 水 上 航行 的 能 力 和 动力 匹配 起 着 
重要 作用 ,直接 决定 了 两 栖 式 清洲 机 器 人 的 工作 效 
率 , 其 核心 问题 是 航行 阻力 的 预测 ” 。 近 年 来 ,不 少 
学 者 通过 计算 流体 力学 ( computational fluid dynam- 
ics ,CFD ) 的 方法 对 浮 性 、 快 速 性 和 其 他 水 动力 特性 
进行 研究 , 周 景 涛 等 ”采用 滑 移 网 格 方法 研究 两 栖 
车 水 上 纵 倾 姿态 的 变化 ,在 该 基础 上 预测 水 上 航行 
阻力 ,并 研究 了 阻力 各 成 分 与 行驶 速度 之 间 的 数学 
关系 。 郑 翔 玉 等 ”通过 静态 网 格 方法 发 现 将 两 栖 车 
辆 的 车 轮 收 起 能 有 效 降低 航行 阻力 。 赵 彬 等 汪 采 用 
混合 而 合算 法 和 动 网 格 技术 研究 静水 直 航 下 两 栖 车 
辆 的 姿态 变化 情况 。 树 宇 峰 等 “通过 重合 网 格 的 方 
法 对 带 有 上 自由 液 面 的 两 栖 车 辆 绕 流 场 进行 分 析 , 并 
得 到 两 栖 车 辆 的 水 动力 学 特性 。 王 少 新 等 ”建立 了 
三 自由 度 的 双 车 厢 两 栖 车 辆 模型 ,并 基于 重奏 网 格 
研究 了 其 航行 阻力 ` 纵 摇 、 垂 荡 运 动 性 能 。 李 海宁 
等 ”采用 重叠 网 格 和 结构 网 格 对 圆柱 的 绕 流 场 进行 
研究 ,分 析 不 同 雷诺 数 下 流体 动力 学 特性 的 多 个 参 
数 的 变化 情况 。 可 见 , 采 用 CFD 数值 计算 方法 研究 
水 动力 特性 已 相对 成 熟 ,本 研究 采用 CFD 数值 计算 
方法 对 两 栖 式 清洲 机 右 人 的 水 动力 特性 进行 研究 。 

本 研究 基于 CFD 以 某 型 两 栖 式 清洲 机 器 人 为 
研究 对 象 ,主要 计算 两 栖 式 清洲 机 器 人 静 浮 态 下 的 


吃水 深度 和 纵 倾角 ,并 预测 静水 直 航 下 的 航行 阻力 。 
首先 通过 重 闪 网 格 方法 对 其 静 浮 态 进行 分 析 , 得 到 
其 静 浮 时 的 吃水 深度 和 纵 倾角 度 ;再 采用 静态 网 格 
和 重 欠 网 格 两 种 不 同方 法 预测 两 栖 式 清洲 机 器 人 在 
静水 直 航 工 况 下 的 航行 阻力 ;最 后 比较 这 两 种 方法 
的 预测 结果 ,分 析 这 两 种 方法 的 适用 性 和 优 劣 性 。 
通过 CFD 数值 计算 得 到 的 两 栖 式 清洲 机 融 人 的 静 
浮 态 和 航行 阻力 相关 数据 为 后 续 机 器 人 的 承载 能 
的 判断 发 动机 的 动力 匹配 、 整 机 的 推进 效率 等 问题 
提供 数据 参考 ,为 掌握 两 栖 式 清洲 机 器 人 的 关键 技 
术 提 供 理 论 依 据 和 技术 支持 。 


1 数学 模型 和 数值 计算 方法 


1.1 研究 对 象 


本 研究 以 某 型 两 栖 式 清洲 机 带 人 为 研究 对 象 ， 
采用 CFD 方法 研究 清 汶 机 器 人 外 流 场 以 及 流 场 和 
清江 机 器 人 之 间 的 相互 作用 。 

简化 的 两 栖 式 清 汶 机 器 人 结构 如 几 1 所 示 , 总 
重 约 为 13.71t, 机 体 长 度 为 10.03 m( 包 含 前 浮 箱 )， 
机 体 宽度 为 5.3 m( 包 含 测 浮 箱 ) , 挖 气 辟 配备 400 L 
铲 斗 ,最 大 航行 速度 为 2 m/s。 


1- 挖 掘 展 ，2- 驾 驶 室 ，3- 发 动机 舱 ，4- 定 位 桩 杜 ; 
5- 机 体 ，6- 侧 浮 箱 ，7- 前 稳定 器 ，8- 前 浮 箱 


图 1 两栖 式 清洲 机 器 人 结构 示意 图 
Fig.1 Structure diagram of amphibious silt cleaner 
基于 CFD 研究 两 栖 式 清 汶 机 器 人 的 水 动力 特 
性 必须 确定 外 流 场 结构 。 保 留 两 栖 式 清洲 机 器 人 最 
外 层 表面 结构 并 生成 面 网 格 , 使 其 成 为 单 连通 区 域 ， 
在 计算 域 和 该 单 连通 区 域 之 间 的 结构 即 为 外 流 场 
结构 。 
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1.2 CFD 数值 计算 方法 


结合 本 研究 使 用 的 CFD 商业 软件 STAR-CCM + 
介绍 相关 数值 计算 方法 。 
1.2.1 控制 方程 和 涡流 模型 

流体 在 流动 过 程 中 遵循 物理 守恒 三 大 定律 : 质 
量 守恒 定律 ` 动 量 守恒 定律 和 能 量 守恒 定律 ,由 于 本 
研究 中 不 涉及 热 交 换 问题 , 故 不 考虑 能 量 守 恒定 律 。 
相应 的 控制 方程 如 下 。 

1) 质量 守恒 方程 :质量 守恒 方程 也 被 称 为 连续 
性 方程 ,其 方程 为 


gp ,akpou) ,dlpv) ,oow) _0 (1) 
ot O% 07 Oz 
式 中 :po 为 流体 的 密度 ;uv、w 为 流体 速度 矢量 在 x、 
yz 轴 上 的 分 量 。 


2) 动量 守恒 方程 :动量 守恒 方程 的 实质 是 牛顿 
第 二 定律 ,本 研究 中 的 流体 均 是 黏 性 为 常数 的 不 可 
压缩 牛顿 流体 , 则 动量 守恒 方程 就 是 Navier-Stokes 
(N-S) 方 程 。 其 公式 为 


工 

p 
式 中 :为 流体 速度 矢量 ;为 哈密 顿 算 子 ; 了 为 单位 
质量 力 ;p 为 压强 ;v 为 运动 黏度 ; ?为 拉 普 拉 斯 
算 子 。 

本 研究 的 汕 流 模型 采用 Realizable k-s 模型 ,该 
模型 既 保留 了 标准 h-s 模型 的 优点 ,如 适用 性 强 , 计 
算 量 适 中 日 精 度 较 好 等 ,又 克服 了 标准 k-s 模型 在 
计算 强 旋 流 和 弯曲 流 时 精度 较 低 的 缺点 ,避免 应 变 
率 较 大 时 计算 出 现 负 应 力 的 问题 ,在 类 似 研 究 的 数 
值 计 算 中 具有 广泛 应 用 。 

1.2.2 计算 域 和 边界 条 件 

计算 域 足够 大 才能 准确 计算 流 场 中 的 物理 量 ， 
但 是 过 大 的 计算 域 会 使 网 格 数量 过 多 ,计算 效率 低 。 
整个 计算 域 边界 如 图 2 所 示 , 入口 设置 为 速度 人 口 ， 
由 清 湾 机 器 人 前 端 向 前 延伸 15 m, 可 在 速度 入 口 设 
置 相 应 的 流体 速度 模拟 航行 速度 ;出 口 设 置 为 静 压 
出 口 ,由 清 淤 机 器 人 后 端 向 后 延伸 30 m; 顶 部 设置 为 
速度 入 口 ,由 清 淤 机 器 人 机 体 中 间 位 置 向 上 延伸 
10 mi; 底部 设置 为 速度 入口 ,由 清 洪 机 器 人 机 体 中 间 
位 置 向 下 延伸 20 m; 左 右 壁 面 设置 为 对 称 面 ,相距 
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域 边界 反射 ,还 需要 将 入 口 .出 口 和 壁面 设置 为 阻尼 
边界 。 


壁面 顶部 


底部 


了 
友 2 ”计算 域 的 边界 条 件 

Fig.2 Boundary conditions of computational domain 
1.2.3 网 格 设置 

本 研究 的 外 流 场 结构 相对 复杂 , 故 采用 局 部 适 
应 性 较 强 的 切割 体 网 格 。 为 了 保证 计算 结果 准确 ， 
在 两 栖 式 清洲 机 器 人 模型 附近 流 场 添加 加 密 区 域 ， 
如 图 3 中 1 区域, 以 获取 清洲 机 器 人 周围 流 场 详细 
的 信息 ;在 水 面 附近 添加 加 密 区 域 ,如 图 3 中 2 区 
域 ,可 获取 准确 的 水 面 波形 。 不 同 区 域 的 网 格 大 小 
不 同 ,为 了 便于 网 格 之 间 的 数据 交换 ,不 同 网 格 之 间 
尽量 采用 2 倍 变 化 率 进行 过 渡 中 。 

为 了 求解 两 栖 式 清洲 机 器 人 在 流 场 中 的 运动 ， 
则 需要 利用 动 网 格 技术 来 满足 运动 过 程 中 网 格 的 变 
化 情况 。 常 用 的 动 网 格 技术 有 网 格 重 构 ` 滑 移 网 格 、 
重 共 网 格 等 ,由 于 重 闪 网 格 技术 适用 于 复杂 模型 的 
网 格 变化 问题 ,因此 本 研究 采用 重 和 至 网 格 技术 , 重 
县 网 格 如 图 3 中 3 区 域 。 


1 区 
FEE 


口 - 帮 
上 加 加 回回 加 加 553533 53 和 
时 上 上 寺村 | Es: 


上 二 


和 rrr i 加 加 面 | 
| 于 上 H 

H Ea 攻占 于 了 HH 
| | 


| | 
[ | | 
HH HE 由 
中 | i 台 加 图 固 名 让 国 | [| | 1 国 | 
上 运 1 本 再 1 | 


3 流 场 网 格 划分 
Fig.3 Mesh generation of flow field 
进行 网 格 无 关 性 验证 后 ,将 本 研究 的 数值 计算 
模型 的 体 网 格 数量 和 质量 整理 如 表 1 所 示 。 
表 1 体 网 格 数量 和 质量 


Tab.1 Quantity and quality of volume mesh 


网 格 数量 0.1 友 网 格 质量 入 17/% 网 格 质量 <0. 1/% 


40 m™ | 。 这 种 设置 方法 能 模拟 出 无 限 宽广 的 流体 空 
间 ,消除 计算 域 边 界 形成 边界 层 造成 的 流 场 素 乱 问 
题 。 为 了 防止 数值 计算 过 程 中 已 产生 的 波 通过 计算 


4 805 752 


1.2.4 自由 液 面 模拟 
自由 液 面 的 模拟 能 提供 航行 过 程 中 水 面 波形 变 


98.885 1.115 
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化 情况 ,也 是 准确 预测 航行 时 兴 波 阻力 的 关键 技术 。 
两 栖 式 清洲 机 右 人 外 流 场 属于 典型 的 气 液 两 相 流 问 
题 ,选择 欧 拉 多 相 分 离 流 物理 模型 定义 水 和 空气 两 
种 流体 。 目 前 模拟 气 液 两 相 流 自由 液 面 的 方法 有 


力 , 静 浮力 的 作用 点 为 排水 体积 的 形 心 , 大 小 为 排水 
的 重量 。 两 栖 式 清洲 机 器 人 水 上 静 浮 问题 从 本 质 上 
来 说 是 重力 和 浮力 的 大 小 以 及 重心 和 浮 心 相对 位 置 


VOF 法 、MAC 法 、Level Set 法 等 ,VOF 法 具有 占用 计 
算 机 资源 少 ,计算 准确 等 优势 ” , 故 采 用 VOF 法 。 

VOF 法 又 称 为 流体 体积 法 ,该 方法 基于 网 格 单 
元 体积 比 的 概念 ,流体 总 体 的 体积 分 数 为 1, 通 过 计 
算 过 程 中 某 种 流体 网 格 所 占 的 体积 比 ,从 而 得 到 各 
项 流体 在 空间 的 分 布 流动 情况 ,进而 确定 自由 液 面 
的 位 置 。 

1.2.5 六 自由 度 体 设 置 

研究 两 栖 式 清 汰 机 器 人 在 水 上 的 运动 情况 时 需 
将 其 设置 为 六 自由 度 体 , 即 治 空间 直角 坐标 系 和 .了 
Z 三 轴 的 平 动 和 绕 三 轴 的 转动 ,从 而 模拟 复杂 的 水 
上 运动 1。 

STAR-CCM + 的 DFBI( dynamic fluid body inter- 
action ) 运 动 模型 实现 了 流 固 动力 学 耦合 ,通过 流体 
和 国体 之 间 的 相互 作用 计算 固体 模型 在 流 场 中 的 运 
动情 况 ,固体 的 运动 亦 影响 到 流体 的 流动 情况 ,如 波 
形 等 。 由 于 两 栖 式 清洲 机 器 人 左右 完全 对 称 , 在 研 
究 浮 性 和 快速 性 时 仅 考虑 升 沉 和 纵 倾 两 个 自由 度 上 


的 运动 。 


2 和 静 浮 态 分 析 


为 了 保证 驾驶 舒适 性 和 视野 开阔 ,通常 要 求 两 
栖 式 清洲 机 器 人 能 水 平静 浮 在 水 面 上 。 通 过 CFD 
的 方法 研究 两 栖 式 清洲 机 器 人 的 静 浮 态 可 得 到 其 准 
确 的 吃水 深度 以 及 纵 倾角 大 小 ,有 利于 调配 两 栖 式 
清洲 机 器 人 的 配 重 ,在 一 定 程度 上 推论 其 承载 能 
计算 得 到 的 静 浮 态 也 是 后 续 各 水 动力 特性 研究 的 
基础 。 


2.1 水 上 静 浮 的 受 力 分 析 


研究 两 栖 式 清洲 机 融 人 更 浮 在 水 面 上 将 建立 以 
下 假设 。 

1) 水 是 不 可 压缩 流体 。 

2) 水 面 平静 、 无 波浪 ,水 面 上 无 风 。 

3 ) 不 考虑 水 面 张力 。 

在 以 上 假设 的 基础 上 ,两 栖 式 清洲 机 絮 人 静 浮 
时 只 受 重力 和 静水 压力 的 作用 。 静 水 压力 在 水 平方 
向 上 的 分 力 相互 抵消 ,垂直 方向 上 的 分 力 则 是 静 浮 


关系 的 问题 ,其 水 上 更 浮 的 平衡 条 件 为 :重心 和 浮 心 
的 作用 线 在 同一 条 垂直 线 上 , 且 重 力 的 大 小 等 于 浮 
力 的 大 小 汪 。 两 栖 式 清 淤 机 器 人 水 上 静 浮 受 力 示意 
图 如 图 4 所 示 ,其 中 0 为 重心 ,6 为 重力 ,C 为 浮 心 ， 
D 为 浮力 。 在 实际 中 ,两 栖 式 清洲 机 带 人 水 上 静 浮 
时 往往 和 水 面 之 间 形 成 某 一 角度 ,该 角度 为 纵 倾角 ， 
即 图 4 中 的 0。 


图 4 两 栖 式 清 淤 机 器 人 水 上 项 浮 受 力 示意 图 


Fig.4 Force diagram of amphibious silt 


cleaner floating on water 
2.2 数值 计算 的 设置 


吃水 深度 是 两 栖 式 清洲 机 器 人 重要 参数 之 一 ， 
由 于 两 栖 式 清 淤 机 器 人 的 形状 并 非 规则 形状 ,所 以 
更 浮 时 的 排水 体积 并 非 关 于 吃水 深度 和 纵 倾 角 0 的 
简单 函数 关系 ,这 就 为 两 栖 式 清洲 机 器 人 水 上 静 浮 
态 的 分 析 和 计算 增加 了 难度 。 

为 了 找 出 两 栖 式 清洲 机 器 人 的 水 上 静 浮 态 ， 
内 相关 研发 单位 常 采 用 水 上 试验 来 确定 。 该 方法 需 
要 生产 出 实物 或 者 比例 模型 ,上 且 需要 配备 相关 精密 
仪器 ,花费 大 且 效 率 低 2” 。 国 内 部 分 学 者 通过 对 
SOLIDWORKS 等 三 维 建 模 软 件 进 行 二 次 开发 ,进而 
研究 两 栖 车 辆 的 静 浮 时 的 位 姿 以 及 相关 的 水 上 性 
能 '"” ,这 些 方法 能 较 准确 地 找到 静 浮 位 姿 , 但 过 程 
繁琐 旦 所 能 求 得 的 水 上 性 能 参数 较 少 ,具有 一 定 的 
局 限 性 。 

本 研究 基于 CFD 研究 两 栖 式 清洲 机 器 人 的 静 
浮 态 ,该 方法 具有 结果 直观 且 花 费 少 等 优势 。 在 上 
文中 介绍 的 CFD 设置 基础 上 ,设置 流速 为 0 的 静 
水 ,两 栖 式 清 汶 机 器 人 水 平 泽 在 水 面 上 (9=0?) , 吃 
水 深度 为 400 mm ,初始 时 刻 两 栖 式 清洲 机 器 人 的 受 
力 情况 如 图 5 所 示 。 
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5 初始 时 刻 的 受 力 图 
Fig.5 Force diagram of initial time 

如 图 5 所 示 : 红 色 矢 量 为 流体 作用 于 两 栖 式 清 
洲 机 器 人 重心 处 的 力 , 其 大 小 等 于 初始 时 刻 的 浮力 
D,; 蓝 色 矢量 为 流体 作用 于 两 栖 式 清洲 机 器 人 重心 
处 的 力矩 。 由 蓝 色 矢量 的 方向 可 知 初始 时 刻 浮 心 C 
的 位 置 在 重心 0 的 前 方 ,此 时 两 栖 式 清洲 机 器 人 受 
力 不 平 衡 ,位 姿 发 生变 化 直到 红色 矢量 和 重力 平衡 、 
蓝 色 矢量 为 0 才 达 到 平衡 状态 。 
2.3 数值 计算 结果 与 分 析 

整个 数值 计算 共和 迭代 6400 步 ,对 应 物理 时 间 约 
为 12.8 s, 记 录 整 个 过 程 中 升 沉 量 和 纵 倾 角 的 变化 
如 图 6 ~7 所 示 。 
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到 6 静 浮 态 下 升 沉 量变 化 曲线 


Fig.6 Curve of heave of center of gravity 
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[7 静 浮 态 下 纵 倾角 变化 曲线 
Fig.7 Curve of pitch angle 
由 图 6 ~7 可 知 ,两 栖 式 清洲 机 器 人 从 初始 位 置 
经 历 了 最 初 的 大 幅度 调整 后 调整 幅度 逐渐 减 小 , 整 
个 过 程 类 似 于 幅 值 减 小 的 三 角 函 数 曲 线 , 并 在 物理 


第 40 卷 


时 间 10 s 之 后 ,两 栖 式 清 淤 机 器 人 基本 达到 了 静 浮 
态 平衡 , 升 沉 量 为 25 mm, 即 静 浮 态 时 吃水 深度 为 
375 mm, 纵 倾角 为 1.6°。 当 液压 臂 前 伸 时 整 机 的 重 
心 前 移 , 浮 心 位 置 也 必须 前 移 以 满足 静 浮 条 件 , 纵 倾 
角 会 减 小 ,因此 保留 少量 的 正 向 纵 倾角 是 很 有 必要 
的 , 即 少量 的 “ 舰 借 " 姿 态 。 同 时 ,少量 的 “ 舰 借 " 密 
态 增 大 了 螺旋 桨 的 浸 深度 ,提高 推进 能 力 。 可 见 ,此 
两 枉 式 清 洪 机 器 人 的 静 浮 态 设计 是 满足 要 求 的 。 

为 了 验证 数值 计算 的 正确 性 ,记录 下 两 栖 式 清 
洪 机 器 人 所 受 浮 力 如 图 8 所 示 。 可 知 浮力 在 大 约 
10 s 后 达到 稳定 状态 , 且 浮 力 最 终 为 134 277. 88 N， 
该 数值 和 两 栖 式 清洲 机 器 人 重力 基本 吻合 。 
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图 8 静 浮 态 下 浮力 变化 曲线 
Fig.8 Curve of buoyancy 


3 航行 阻力 预测 


两 栖 式 清洲 机 器 人 水 上 航行 阻力 对 其 动力 匹配 
有 着 十 分 重要 的 指导 作用 。 此 外 ,合理 分 析 清 洲 机 
器 人 的 航行 阻力 有 助 于 对 其 进行 减 阻 优化 ,提高 其 
航行 速度 ,加 快 工作 效率 ,同时 有 效 控制 污染 排放 及 


3.1 航行 阻力 成 分 


当 两 栖 式 清洲 机 器 人 在 水 上 匀速 航行 时 ,其 航 

行 阻力 分 为 水 阻力 和 空气 阻力 ,一般 认 为 空气 阻力 

太 小 而 忽略 。 水 阻力 及 是 摩 所 阻力 R;、 形 状 阻力 

Rw 和 兴 波 阻力 Rw 之 和 ,其 中 形状 阻力 和 兴 波 阻力 
也 被 称 为 剩余 阻力 R; ,其 关系 式 为 

Ri = Re + Roy + Ry (3) 

R= Roy + Rw (4) 

式 中 :摩擦 阻力 是 机 体 与 水 之 间 的 摩擦 力 , 即 运动 方 

向 上 的 切 应 力 ; 形 状 阻力 是 清洲 机 絮 人 前 后 表面 压 

力 不 等 所 形成 的 阻力 ,其 成 因 是 流体 沿 着 机 器 表面 

流动 时 发 生 在 边界 层 的 分 离 现 象 ; 兴 波 阻力 是 由 于 

波浪 导致 压力 分 布 改 变 而 产生 的 阻力 ;剩余 阻力 均 

是 由 压 差 导致 的 ”。 
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3.2 基于 静态 网 格 的 航行 阻力 预测 


在 低速 航行 时 航行 姿态 对 于 航行 阻力 影响 较 
小 ' ,因此 考虑 基于 静 浮 态 采用 静态 网 格 来 预测 航 
行 阻力 。 

根据 静 浮 态 分 析 中 得 到 的 数据 对 模型 位 姿 进 行 
调整 ,采用 静态 网 格 预测 两 栖 式 清洲 机 如 人 最 高 航 
速 (2 m/s) 时 静水 直 航 工 况 下 的 航行 阻力 。 当 预测 
的 航行 阻力 数值 趋 于 稳定 后 停止 数值 计算 ,整个 数 
值 计算 过 程 共 迄 代 3 200 步 , 该 过 程 中 的 阻力 变化 如 
图 9 所 示 。 


航行 阻力 / 


0 400 800 1200 1600 2000 2400 2800 3200 


过 代 数 
图 9 静态 网 格 下 航行 阻力 变化 曲线 
Fig.9 Curve of sailing resistance 
由 图 9 可 知 ,航行 阻力 在 迭代 1 400 步 后 趋 于 稳 
定 ,最 终 数值 为 3405.58 N ,其 中 水 阻力 为 3396.77 NN， 
空气 阻力 为 8. 83 N, 由 于 空气 阻力 过 小 , 故 不 予 考 
虑 。 水 阻力 系数 C, 的 计算 公式 为 
Fi 
CT (5) 
式 中 :下 为 水 阻力 ;p 为 水 的 密度 ;v 为 航速 ;5 为 两 
栖 式 清 淤 机 器 人 水 下 的 湿 表 面积 ” 。 计 算 两 栖 式 清 
湾 机 器 人 的 水 阻力 系数 C, 为 0.029。 各 部 件 的 水 阻 
力 成 分 如 表 2 所 示 。 
表 2 ”静态 网 格 下 各 部 件 的 水 阻力 


Tab.2 Water resistance of different parts 


部 件 摩擦 阻力 AN 剩余 阻力 AN 合计 /N 
机 体 180. 65 1 406.72 1 587. 37 
前 稳定 器 44. 94 1 764. 46 1 809.40 
合计 225,59 3171.18 3 396.77 


从 水 阻力 成 分 看 ,摩擦 阻力 占 航 行 总 阻力 的 
6. 64% ,剩余 阻力 所 占 比例 为 93.36% ;从 各 部 件 所 
受 水 阻力 大 小 看 ,前 稳定 融 所 受 水 阻力 占 比 达到 
53.27% ,整个 机 体 ( 包 括 前 浮 箱 和 测 浮 箱 ) 所 受 的 
水 阻力 占 比 为 46.73% 。 


式 清 次 机 器 人 水 动力 特性 研究 
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3.3 基于 重 人 释 网 格 的 航行 阻力 预测 


3.3.1 航行 阻力 预测 结果 

采用 重 盖 网 格 计算 两 栖 式 清洲 机 妖 人 2 m/s 航 
速 下 由 设置 的 初始 状态 达到 平衡 状态 过 程 的 升 沉 
量 、 纵 倾角 以 及 航行 阻力 的 变化 情况 。 整 个 数值 计 
算 共 和 迭代 15 000 步 ,对 应 物理 时 间 为 30 s, 整 个 过 程 
中 航行 阻力 如 图 10 所 示 。 


0 2 4 IO 16 18 20 9224 26 28 30 
物理 时 间 /s 


图 10 重合 网 格 下 航行 阻力 变化 曲线 
Fig. 10 Curve of sailing resistance 
由 图 10 可 知 ,航行 阻力 于 10 s 后 趋 于 稳定 ,最 
终 数 值 为 3594.40 N ,其 中 水 阻力 为 3585.84 N , 空 
气 阻力 为 8.56 N, 由 式 (5) 计 算得 到 两 栖 式 清洲 机 
器 人 的 水 阻力 系数 Cu 为 0.031。 不 同 部 件 所 受 的 水 
阻力 如 表 3 所 示 。 
表 3 重生 网 格 下 各 部 件 的 水 阻力 


Tab.3 Water resistance of different parts 


部 件 摩擦 阻力 /N 剩余 阻力 /N 合计 /N 
机 体 191.14 1284.41 1 475. 55 
前 稳定 器 44. 60 2 065. 69 2 110.29 
合计 235.74 3 350. 10 3 585. 84 


从 水 阻力 成 分 角度 看 ,摩擦 阻力 占 水 阻力 的 
6. 57% ,剩余 阻力 占 比 为 93. 43% 。 该 结论 与 文献 
[14] 的 结论 “在 伟 如 德 系数 Fy 较 小 阶段 (Ff,。 < 
0. 30) ,两 栖 车 航行 阻力 以 黏 压 阻 力 为 主 , 占 比 为 
60% ~70% ,摩擦 阻力 在 10% 左右 , 兴 波 阻力 通常 
不 超过 30% ” 妆 基 本 一 致 ,这 是 由 于 两 栖 式 清 淤 机 
器 人 外 形 为 箱 型 体 ,在 航行 时 会 产生 大 量 的 流动 分 
离 和 兴 波 ,导致 剩余 阻力 很 大 ,而 低速 航行 时 摩擦 阻 
力 较 小 。 

从 各 部 件 所 受 水 阻力 大 小 看 ,前 稳定 器 的 水 阻 
力 占 比 为 58. 85% , 可见 前 稳定 器 对 于 两 栖 式 清 淤 
机 器 人 的 航行 阻力 影响 较 大 。 

以 重 公 网 格 数值 计算 的 结果 为 基础 ,计算 静态 
网 格 方法 所 得 各 项 阻力 的 偏差 ,统计 如 表 4 所 示 。 

由 表 4 可 知 ,以 静 浮 态 为 基础 采用 静态 网 格 对 
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两 栖 式 清洲 机 右 人 航行 阻力 的 预测 结果 和 采用 重 伙 
网 格 方法 的 预测 结果 整体 相差 较 小 ,但 对 于 前 稳定 
器 的 剩余 阻力 预测 结果 相差 较 大 。 

表 4 各 阻力 的 偏差 


Tab.4 Deviation of resistance 


阻力 成 分 遍 差 /( % ) 
机 体 摩擦 阻力 -5.49 7.58 
剩余 阻力 9.52 
前 稳定 器 摩擦 阻力 0.76 -14.26 
剩余 阻力 -14.58 
整体 摩擦 阻力 -4.31 -5.27 
剩余 阻力 -5.34 
空气 阻力 3.15 
3.3.2 浮力 及 位 姿 变 化 情况 


记录 下 整个 过 程 中 浮力 的 变化 情况 如 图 11 
所 示 。 浮 力 于 10s 后 趋 于 稳定 , 最 终 数值 为 
134 208. 86 N ,该 数值 和 两 栖 式 清洲 机 器 人 重力 基本 
钨 合 。 
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图 11 航行 态 下 浮力 变化 曲线 

Fig. 11 

由 于 两 栖 式 清洲 机 器 人 水 上 航行 是 一 个 动态 过 

程 ,因此 其 位 姿 变 化 处 于 动 平衡 。 整 个 过 程 两 栖 式 

清洲 机 带 人 的 升 沉 量 和 纵 倾角 变化 情况 如 图 12 ~ 
13 所 示 。 


Curve of buoyancy 
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图 12 航行 态 下 升 沉 量变 化 曲线 
Fig.12 Curve of heave of center of gravity 
可 图 12 ~13 可 知 , 航 行 过 程 中 的 升 沉 量 和 纵 倾 
角 变 化 较 小 ,这 是 由 于 两 栖 式 清洲 机 器 人 重量 大 且 
航速 低 。20 s 之 后 清 淤 机 器 人 基本 处 于 动 平衡 的 状 
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态 , 升 沉 量 为 -4 mm, 即 吃水 深度 为 404 mm , 纵 倾角 
为 1.19*。 发 现 两 栖 式 清洲 机 右 人 低速 航行 时 会 出 
现 吃 水 深度 增加 、 纵 倾角 减 小 , 即 机 融 下 沉 、 出 倾 程 
度 增 大 的 规律 。 这 是 由 于 水 流 流 过 机 器 时 使 得 流速 
增加 ,从 而 导致 水 作用 于 机 顺 的 压力 减 小 ,机 融 略 微 
下 沉 , 且 机 首 处 水 流速 度 增 大 更 多 ,导致 机 器 出 倾 程 
度 增加 , 机体 底 部 的 压力 云图 (图 14) 可 验证 该 
规律 。 
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图 13 航行 态 下 纵 倾 角 变 化 曲线 
Fig. 13 Curve of pitch angle 
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图 14 机 体 底部 压力 云图 
Fig. 14 Pressure cloud chart of bottom of machine body 

3.3.3 前 稳定 器 的 水 阻力 分 析 

前 稳定 器 相对 于 清 淤 机 器 人 机 体 而 言 体 积 较 
小 ,但 受到 了 较 大 的 水 阻力 , 且 采 用 静态 网 格 和 重 县 
网 格 两 种 方法 对 其 剩余 阻力 预测 结果 相差 较 大 , 现 
对 前 稳定 带 的 水 阻力 进行 分 析 。 

根据 前 稳定 器 前 后 压强 图 (图 15 ) 发 现 :两 栖 式 
清 淤 机 器 人 的 前 稳定 器 在 其 A、B 区 域 处 前 后 高 压 
面积 和 压强 相差 较 大 ,这 直接 导致 了 前 稳定 器 受到 
较 大 剩余 阻力 。 由 波形 高 度 图 (图 16) 可 知 :A 区 域 
主要 是 由 于 流 过 前 浮 箱 的 高 速 流体 冲击 造成 前 后 压 
差 较 大 ;B 区 域 则 是 由 于 在 该 处 产生 了 类 似 于 两 栖 
车 辆 车 首 波 系 的 波浪 ,导致 此 处 高 压 面积 和 前 后 压 
差 相 差 较 大 。 综 上 所 述 ,前 稳定 器 受到 的 剩余 阻力 
较 大 。 又 由 于 前 稳定 器 的 形状 圆滑 无 机 角 , 所 以 水 
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Fig. 1$ Pressure diagram of front stabilizer 
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图 16 航行 时 波形 高 度 图 
Fig.16 Wave height while sailing 
造成 前 稳定 器 剩余 阻力 较 大 的 原因 复杂 , 且 吃 
水 深度 和 纵 倾 角 能 直接 影响 到 高 压 面 积 的 分 布 情 
况 ,静态 网 格 方法 基于 融 浮 态 而 忽略 了 航行 时 两 栖 
式 清洲 机 絮 人 的 位 姿 变 化 ,因此 会 造成 较 大 的 偏差 。 


4 结 论 


本 研究 采用 重大 网 格 对 某 两 栖 式 清洲 机 器 人 水 
上 静 浮 态 进行 数值 计算 并 得 到 其 吃水 深度 和 纵 倾 
角 ; 采 用 不 同 的 方法 对 两 栖 式 清 淤 机 器 人 水 上 航行 
阻力 进行 预测 ,对 比 并 分 析 了 两 种 方法 的 优 劣 。 本 
研究 的 计算 结果 可 为 后 续 水 中 推进 .发 动机 的 动力 
匹配 等 问题 提供 数据 支持 ,为 掌握 两 栖 式 清 淤 机 器 
人 的 关键 技术 提供 理论 依据 和 技术 支持 。 得 到 以 下 
结论 。 

1) 通过 数值 计算 得 到 某 型 两 栖 式 清洲 机 器 人 吏 
浮 态 下 的 吃水 深度 为 370 mm, 纵 倾角 为 1. 6"; 航 速 
为 2 m/s 时 静水 直 航 的 航行 阻力 预测 为 3 594. 4 N， 
水 阻力 系数 为 0.036。 

2) 采 用 静态 网 格 和 重 秋 网 格 均 能 预测 两 栖 式 请 
与 机 器 人 航行 的 航行 阻力 。 从 计算 效率 看 ,前 者 大 
代 3 200 步 即 可 得 到 稳定 的 计算 结果 ,计算 效率 较 
高 ;后 者 迭代 15 000 步 可 得 到 稳定 的 计算 结果 ,计算 
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效率 较 低 。 从 适用 性 看 ,前 者 由 于 忽略 了 两 栖 式 清 
六 机 器 人 水 上 航行 姿态 的 变化 , 仅 适 用 于 航行 姿态 
变化 较 小 的 低速 航行 工 况 ;后 者 能 计算 航行 姿态 的 
变化 情况 ,适用 于 各 种 航行 工 况 。 

3 ) 两 栖 式 清 汶 机 器 人 低速 航行 时 会 出 现 机 器 下 
沉 、 出 倾 程度 增 大 的 规律 ,这 主要 是 由 于 机 器 附近 水 
流 流速 增加 ,水 作用 于 机 咒 的 压力 减 小 造成 的 。 
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